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Abstract 

We complete the next-to-next-to-leading order QCD calculation of the branching ratio 
for B — > X s £ + £~ including recent results for the three-loop anomalous dimension matrix 
and two-loop matrix elements. These new contributions modify the branching ratio in the 
low-g 2 region, BR«, by about +1% and —4%, respectively. We furthermore discuss the 
appropriate normalization of the electromagnetic coupling a and calculate the dominant 
higher order electroweak effects, showing that, due to accidental cancellations, they change 
BRu by only —1.5% if a(fi) is normalized at fi — 0(m b ), while they shift it by about —8.5% 
if one uses a high scale normalization \x = 0(M W ). The position of the zero of the forward- 
backward asymmetry, q$, is changed by around +2%. After introducing a few additional 
improvements in order to reduce the theoretical error, we perform a comprehensive study 
of the uncertainty. We obtain BR«(1 GeV 2 < q 2 < 6GeV 2 ) = (1.57 ± 0.16) x 10" 6 and 
g 2 = (3.76 ± 0.33) GeV 2 and note that the part of the uncertainty due to the fe-quark mass 
can be easily reduced. 



1 Introduction 



The rare semileptonic b — > s£ + £~ transitions have been recently observed for the first 
time by Belle and BaBar in form of the exclusive B — > K^*'£ + £~ modes [1], and also 
inclusive measurements are now available [2] . Like all other flavor-changing- neutral-current 
transitions, these channels are important probes of short-distance physics. Their study 
can yield useful complementary information, when confronted with the less rare b —* S7 
decays, in testing the flavor sector of the Standard Model (SM). In particular, a precise 
measurement of the inclusive B — > X s £ + £~ channel would be welcome in view of new 
physics searches, because it is amenable to a clean theoretical description, especially for 
dilepton invariant masses, m 2 u = q 2 , below the charm resonances, namely in the range 
1 GeV 2 < m 2 u < 6 GeV 2 . In theoretical calculations the reference window is usually 0.05 < 
s = m 2 u /ml < 0.25. 

The SM calculation of the Branching Ratio for B — > X s £ + £~ in this low-s region, 
BR« in the following, has reached a high degree of maturity [3-6]: it now includes a 
nearly complete resummation of large QCD logarithms L = In mb/M w up to Next-to- 
Next-to-Leading Order (NNLO). The peculiarity of the b — > s£ + £~ amplitude is that it 
involves large logarithms even in the absence of QCD interactions. Factoring out G^a, 
it receives contributions of 0(a™L n+l ) at Leading Order (LO), of 0(a™L n ) at Next-to- 
Leading Order (NLO), and of O^ot^U 1 " 1 ) at NNLO. Since the QCD penguin operators 
Qs-Qe are suppressed by small Wilson coefficients, which makes the computation of their 
two- loop matrix elements dispensable 1 , the only potentially relevant NNLO terms still 
missing at low s originate from the three-loop Anomalous Dimension Matrix (ADM) of the 
operators in the low-energy effective Hamiltonian, and from the two-loop matrix element 
of one of them, namely the vector-like semileptonic operator Q 9 . Finally we would also 
like to remark that the NNLO matrix elements of the current-current operators Q± and Q2 
for the interesting high-s tail of the B — > X s £ + £~ decay distributions have been calculated 
just recently [8]. 

On the other hand, electroweak effects in b — > s£ + £~ have never been studied in the 
literature. As shown in the case of radiative decays [9-11], they may be as important as the 
higher order QCD effects. BR« is generally parameterized in terms of the electromagnetic 
coupling a, but the scale at which a should be evaluated is, in principle, undetermined 
until higher order electroweak effects are taken into account. This has led most authors 
to use a{mu) « a^rrif,) w 1/133. Indeed, as we will discuss later on, in the absence of 
an 0(a) calculation, there is no reason to consider a(m«) more appropriate than, say, 
a(M w ) ~ 1/128. As BR^ is proportional to a 2 , the ensuing uncertainty of almost 8% is 
not at all negligible. 

The purpose of this paper is threefold: i) in Section |2 we evaluate the two main missing 
NNLO QCD contributions to B — > X s £ + £~ taking advantage of the new results for the 
three-loop ADM [12,13] and of existing calculations of 0(a 2 ) corrections to semileptonic 

1 ln B — > X s -f the two-loop matrix elements of Q3-Q6 reduce the branching ratio by around 1% [7]. 
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quark decays; ii) in Section |3] we study to what extent higher order electroweak effects 
influence this decay mode by calculating the dominant 0(a) contributions to the running 
of the Wilson coefficients and estimating other potentially large effects; Hi) in SectionEJwe 
update the SM prediction of BR« and of the position of the zero of the Forward-Backward 
(FB) asymmetry, sq, using the above results, and discuss a few improvements aimed at 
reducing the residual theoretical error in those observables. 

2 Completing the NNLO QCD Calculation 

Let us first focus on the missing NNLO QCD contributions. The main steps of the NNLO 
calculation have been outlined in [3] , while two-loop matrix elements and bremsstrahlungs 
corrections have been computed in [4] . In the following we adopt the operator basis of [3] , 
which we enlarge to include the electroweak penguin operators Q^-Qq and the operator 
Q\, in view of our subsequent discussion of electroweak effects: 



Qi = 


(s Ll ,T a c L )(c L ^T a b L ) , 


Q 2 = 


0l7m c £)( c £7^l) , 


Q 3 = 


(sL7A)E g (<?7 M <?) , 


Q 4 = 


(s L l,T a b L )j: q (qYT a q) , 


Q 5 = 






(sLl»lvl P T a b L )Z g (qYl u YT a q) , 


Q? = 


(s L%l b L )J2 q Q q (qYq) , 


Q 4 Q = 


{Hl»T a b L )Y. q Q q {<tfT a q) , 


Q Q ,= 


(s L %^%b L )E q Q q (qYi u Yq) , Qe = 


(s L i^ P T% L )T, q Q q (qYYYT a q) , 


Q\ = 




Qi = 


^m b (s L a^b R )F^, 




j-m b (s L a^T a b R )G« u , 


0,9 = 


S(*L7A)£#7^) , 


Qio = 


^(sLl,b L )Z e (hV) ■ 



The operator Q\ corresponds in four dimensions to (sl7 m 6l)(6i,7 m 6l) and receives contri- 
butions from electroweak boxes [14] . In order to avoid the introduction of traces with 75 
at all orders in QCD, we have rewritten this operator in a slightly more involved form 2 . 

The renormalization scale dependence of the Wilson coefficients C T (fi) = (Ci(/i), . . . , 
Cio(/f)) of the effective operators is governed by the Renormalization Group Equation 
(RGE) whose solution is schematically given by 

C(fj) = U(fi,iio,oi)C(fio) ■ (2) 

In the case at hand we are interested in the running of the Wilson coefficients from the 
electroweak scale /xo = 0(M W ) down to the low-energy scale fib = 0(m&). Neglecting 
for the moment the electromagnetic coupling a, all the QCD corrections to the initial 
Wilson coefficients relevant at NNLO have explicitly been given in [3], while the only 

2 We thank M. Misiak for this suggestion. 
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di 
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0.0107 
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0.0071 
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-0.0008 
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fi 


0.0010 


0.0130 


0.0045 


-0.0022 


-0.0714 


-0.0008 


0.0299 
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0.0035 








hi 


0.0114 


-0.0107 


-0.0012 


-0.0057 


-0.0098 


0.0002 


0.0122 









Table 1: Numerical coefficients parameterizing the RGE solutions in Eqs. (J3J), (ITU) and (IT2"|) . 

0(a 2 ) contributions to the evolution matrix U(fi, Ho, a ) pertinent to b — > s£ + £~ at NNLO 
concerns the mixing of Q2 into Q-j and Qg. Expanding U(fj,,/j, Q , a) in powers of « s (/ifc)/47r 
- see Eq. © — we denote these terms by £vm/-*b, Ho) an d U^g 2 (fib, Ho), respectively 3 . The 
ingredients necessary for the calculation of Uf^ 2 Ho) were already available four years 

ago, and are included in [3]. On the other hand, the calculation of £^92 (/•*&, Ho) requires 
the knowledge of the three-loop self-mixing of Q\-Qe as well as the three-loop mixing of 
Qi-Qg into Qg. The relevant ADM entries have been calculated only recently in [12,13]. 
Solving the NNLO RGE [13, 15], we obtain 

U%(Hb: Ho) = £ [hv ai + Ci V ai+1 + d lV a > +2 ] « 181.42 e^ 1901 " r? 1 ' 5212 , (3) 

i=l 

where 77 = a s (/xo) / a s (/Xb), and the so-called magic numbers Oj, bi, Ci and di can be found in 
Table [TJ Using a s (M z ) = 0.119 and mt, = 4.8 GeV, we find numerically Ug Q2 (mb,M w ) ps 
4.1. Unless 77 < 0.48, our result is within the range that was guessed in [3], —10 77 < 
Us miHb, Ho) < 10 r]. Our determination of 92 (/^b, /^o) eliminates one source of uncertainty 
in the NNLO calculation — we will study the residual scale dependence later on — and 
increases the value of BR« presented in Section|3]by about 1%, the exact amount depending 
on the choice of the various renormalization scales. 

Another missing ingredient of the complete NNLO calculation is the two-loop 0(a 2 ) 
matrix element of Qg. This contribution is necessary because Qg has a non- vanishing 
matrix element at tree-level as well as a non- zero Wilson coefficient in leading logarithmic 
approximation. Fortunately, no explicit calculation is necessary in this the QCD 

corrections to b — > s£ + £~ are identical to those to b — ► u£v (t — * b£u), in the limit of 
vanishing strange (bottom) quark mass. In particular, we want to have the QCD corrections 

3 The definition of Ug^ (fib, Mo) adopted in Eq. (37) of Rcf . [3] corresponds to rjUg 92(^6; Mo) used here. 
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Figure 1: Second order perturbative corrections to the matrix element of Q g , see Eq. flU), as 
a function of the dilepton invariant mass. The lower (upper) red curve corresponds to the 
1 — s expansion up to third (fourth) order [16], while the blue curves correspond to the s 
expansion [17] (central value and linearly added errors). The shaded areas show the cuts of the 
experimental analysis of the Belle and BaBar Collaborations [2] to reject backgrounds due to 
J/vp resonances. 

to the dilepton invariant mass spectrum. The 0(a 2 ) corrections to this spectrum for the 
decay b — > uiv have been computed in [16] in terms of an expansion in 1 — s, which 
converges also well in the low-s region of interest. The results up to third and fourth order 
in the 1 — s expansion are shown in Figure EJ in the 6-quark pole mass scheme. 

The dilepton invariant mass spectrum at small s can also be obtained from the M^/M 2 
expansion of the second order QCD corrections to the top quark decay calculated in [17]. 
In principle, using the replacements M w — > m« and M t — > nib, this expansion is better 
suited for the low-s region. However, Ref. [17] provides only the terms up to M^/M/, 
obtained through a Pade approximation from a q 2 /M 2 expansion. The ensuing uncertainty 
is displayed in Figure^ where the errors of the various coefficients have been added linearly 
and as before the pole mass scheme for the 6-quark mass has been used. Although this is 
likely to overestimate the uncertainty in the calculation based on [17], the precision achieved 
is sufficient for our purposes. Moreover, Figure [T] shows that the two approaches [16,17] 
to the calculation of the 0(a 2 ) corrections to the invariant dilepton mass spectrum agree 
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Figure 2: Some of the basic diagrams contributing to b — > s£ + £ in the effective theory. 



quite well in the low-s region. Normalizing the matrix element to its tree level value, 

2 



(Q9) 



(0) I 2 



1 



a s (m b ) (i)^^ t ( a s (m b ) \ (2 ) 



7T 



-Ogg (S] 



(4) 



a simple approximation of the real and virtual corrections denoted by a gg (s) is 

4 2 9 } (s) « -19.2 + 6.1 s+ (37.9 + 17.21ns) s 2 - 18.7 s 3 . (5) 

This approximate formula is valid in the range < s < 0.4, in the pole mass scheme for the 
6-quark mass. Using this expression in the NNLO calculation of BR« presented in Section 
EJ we observe a reduction of about 4%, that overcompensates the three-loop running of 
Cg(fib) discussed above. 



3 Higher Order Electroweak Effects 

A comprehensive study of electroweak effects has been performed in the case of B — > X s ^ 
[9-11] but, to the best of our knowledge, such corrections have never been considered in 
the calculation of BR«. Unlike in the b — > 57 processes, only virtual photons appear in 
the basic b — > si + l~ transition. In the following, we will try to identify and estimate the 
dominant 0(a) electroweak effects. When we consider large logarithms we refer only to 
terms like hnm b /M w and exclude the case of very small s, which would also involve large 
logarithms of the form In s. 

Let us first concentrate on the origin of large logarithms at 0(a). To this end, it is 
best to work from the very beginning in the effective theory at lowest order in the electro- 
magnetic coupling. At the weak scale, we distinguish between two kinds of contributions 
to the b — > s£ + £~ amplitude: a) those containing a virtual photon and b) those that do 
not include photons at all. Examples of the first kind are Q2 and Q? insertions, while Qg 
and Qio insertions belong to the second category — see Figure El This classification holds 
true also at higher orders in QCD. 

Evolving from higher to lower scales, the contributions of class a) receive 0(a) correc- 
tions involving large logarithms related either to the running of the QED coupling constant 
or to the evolution of the relevant Wilson coefficients. The first effect has its origin in vac- 
uum polarization diagrams on the virtual photon propagator proportional to alnm«//i, 
like the one shown in Figure E3 that can be effectively taken into account by changing the 
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Figure 3: Vacuum polarization diagrams. 



scale at which a is evaluated in the overall factor of BR« from the weak scale to ma- The 
second effect has to do with the fact that, in full analogy to what happens in QCD, pho- 
tonic interactions induce an ADM for the effective operators — see Figure EJ In order to 
include all 0(a) logarithms, it is therefore sufficient to evolve both the photonic coupling 
and the Wilson coefficients according to a mixed QED-QCD RGE from the matching scale 
/i = 0(M W ) to a low-energy scale \i = 0(m b ). Beyond leading order in a, one can simply 
normalize the amplitudes of type a) with G^a(fi), and evolve the Wilson coefficients at a 
fixed, as done in [9, 11]. 

Concerning contributions of class b), no vacuum polarization diagram can contribute at 
0(a), and the only 0(a) large logarithms are related to the running of the operators. As 
the normalization of the coupling constant is determined by the short distance interactions 
at the matching scale, the Wilson coefficients are appropriately expressed in terms of 
using the relation g 2 (M w )/4ir = a(M w )/ sin 2 9 W (M W ) = V2/n G^M^, where all running 
couplings are in the MS scheme. The typical case is that of box diagrams, which are clearly 
proportional to g 4 (M w ) / oc G 2 M^. After decoupling, all short distance information is 
encoded in the Wilson coefficients and in G^, which does not evolve in the effective theory. 
The coupling e 2 in front of Qg and Qio — see Eq. ([Q — is naturally evaluated at the weak 
scale. 

All electroweak couplings entering the matching can be expressed in terms of G^, M w , 
and sin^^M,^). These parameters are frozen in the low-energy effective theory. On the 
other hand, the electromagnetic coupling of the photon in the low-energy effective theory 
does run and the low-energy matrix elements must be expressed in terms of the low-energy 
gauge couplings a s (n) and a([i), and of the Wilson coefficients C(jj). This simple intuitive 
picture is naturally realized by a mixed QED-QCD RGE. In fact, a general feature of 
the Renormalization Group (RG) evolution is to keep the coupling in front of the Wilson 
coefficients fixed at the high scale via factors of a(fj, )/a(/u l ). 

The only complication we face here is due to our choice of the operator basis. The factor 
e 2 / g 2 in front of the operators Qg and Qio is such that the Wilson coefficients at the high 
scale Cg(^o) and Cio(fio) start at 0(a s G fl ) rather than at 0(aG Pi ). As a consequence, the 
evolution of the coefficients Cg(/i) and Cio(fi) does not explicitly include a rescaling factor 
a(fio) / a([i) as, for instance, in the standard QCD evolution of Ci(ji) and C 2 (/x). However, 
we will show at the end of this section that the correct normalization of the electroweak 
coupling in front of these contributions is preserved by the QED (3 function coefficients 
(3^ and that appear in the self-mixing of Qg and Qio — see Eqs. (J3~%]) and (|3~9"j) — due 
to the e 2 factor in the definition of these operators. 
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Figure 4: Some diagrams contributing to the running of the Wilson coefficients. 

Summarizing, we perform the RG evolution of the Wilson coefficients which results 
in the normalization of the whole amplitude in terms of a(fi). As we will discuss later 
on, additional sub-dominant corrections come from 0(a) matrix elements of the effective 
operators and from 0(a) corrections to the matching conditions. We employ sin 2 6 W = 
sm_ 2 ^f (M z ) = 0.2312, M t = M t (M t ) = 165 GeV and M w = 80.426 GeV. We adopt an 
MS definition of a(ji) with a {5 \M z ) = 1/127.765 [19] as input. 

Let us now study explicitly the QED-QCD evolution of the Wilson coefficients. Our 
basis (J]) differs from the one used in [9] by a factor g~ 2 in the normalization of Q7- 
Qio, which complicates slightly the counting of couplings and the comparison with the 
latter paper. Neglecting the running of a — an 0(a 2 ) effect in the running of the Wilson 
coefficients — we can expand the evolution matrix U /i , a) of Eq. (J2J) in the following 
way 



%Po,a) = t>i° Wo) + ^f^Wo) + (^)VWo) 

+ ^rt^Wo) + f t/ e (1 Wo) + ^^ e {2 Wo) + • • • • 
a s (n) 47r {4ir) z 



(6) 



The matrices f/i fio) an d Ue (/■*> A*o) describe pure QCD and mixed QED-QCD evolu- 
tion, respectively, and are functions of the ADM 4 



^(0) , 0^ ,(1) 

4tt e (4tt) 2 Je 



(7) 



of the operators in question and of the QCD, QED and mixed QED-QCD (3 functions. 
The matrices (/i, /io) and U^ifi, /i ) can be computed using the formalism described 
for instance in [20]. Expanding also the Wilson coefficients at the weak scale 



CM = c?Vo) + ^c?«M + (^) 2 ^ 2) (/^o) 

+ ^ 1)( ^ )+ w^ 2)M+ ---' 



(8) 



4 Eq. J7J) contains also a term proportional to a 2 /a s (/i), from the mixing of Qg into Q3 . We have 
checked that its contribution to the evolution of the relevant Wilson coefficients is negligibly small. 
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and inserting © and JSJ) into ©, we find the expressions for the various terms at the low 
scale /i 

Here the coefficients Oi (/j) result from the 0(1), 0(a s ) and 0(a 2 ) contributions to C(/iq) 
and from the QCD evolution matrices £/i^ (/i, /io), whereas the Ci^ (/i) stem from the various 
0(a) corrections to O(/z ) and the QCD, QED-QCD evolution matrices Z/i^ (/x, /x ) and 
£> e W (/i,/i ). 

The formally leading electroweak effects are 0(aa™ _1 L n+1 ) in the amplitude and are 
contained in 

C?^) = ^e (0 Wo)C?^o) • (10) 

The non- vanishing 0(a) mixing described by 7e°^ has been calculated in [10], except for the 
QED mixing of Qg and Qio- The lowest order — 0(a s ) in our notation — mixing between 
the electroweak penguin operators Q®-Qq and Qg is also missing in the literature, as is 
the 0(a s ) mixing of the operator Q\. We have computed all missing entries and report 
the full matrices 7e°' ) and in the Appendix. Since only the mixing of Q2 into Qg and 
Qio is relevant at this order, we solve the RGE and get 

d%) = £ [e^- 1 + fa*] « 5.7947 e" 11 - 508 ^ 2 - 8808 , (11) 

i=l 

where we used (/-to) = 1, and the magic numbers Oj, ej and /j are given in Table HJ The 
approximation is valid within less than 1% for 0.5 < r] < 0.6. The formally leading QED 
contribution shifts Cg{m h ) by only 0.00006 for 77 = 0.56. Since Qg mixes into Qio at O(a) 
we also have a contribution to 

9 

= + M*] « 1-9013 e"™. ^1.2089 ? (12) 

t=l 

which shifts Cio(77i&) by 0.00014 for 77 = 0.56. Again the numerical coefficients a^, Qi 
and /tj can be found in Table d and the approximation is valid within less than 1% for 
0.5 <r]< 0.6. Since in our normalization Og(m&) ~ 0.072 and Oio(m&) ~ —0.072 at NNLO 
QCD, the impact of these 0(a/a s ) corrections on BR^ is tiny. 

As discussed in [9,21], next-to-leading electroweak effects can be larger than the leading 
ones. Moreover, in the case of b — > s£ + £~ , the LO QCD contribution is accidentally small 
compared to the NLO QCD contribution. Let us therefore investigate the importance of 
these 0(aa™ _1 L n ) effects. The general expression for Ce(fj) is 

0? } (/i) = fW(jM,i*)dV>Qio) +rjUi 0) (»^o)Ci 1 \» ) + U?\^^)C?\^) - (13) 
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The last term in this equation requires the knowledge of the 0(aa s ) ADM ^\ which we 
have calculated using the method described in [12]. The result is given in the Appendix. 
Rather than giving analytic expressions for these 0(a) corrections to the Wilson coefficients 
of Qr, Q9 and Qio, we only provide some approximate formulas and refer to a forthcoming 
publication [22] for more details. Taking into account the different operator normalization, 
we find the same results as in [9] for the 0(a) correction to the Wilson coefficient of the 
electromagnetic operator, 

Cf 7 V) « 0.2894 e -5.1605, 1.4543 n.3060 + 1Q654 ^14.857, 8.2193 

« 0.0107 + 0.0221, 

where we have distinguished between the last two terms on the right-hand side of Eq. ()13|) 
- in the particular case of CV(/i) the first term is exactly zero as there is no 0(a s ) 
mixing of the four-quark operators into the magnetic one, and the 0(a) contribution to 
Cr(no) also vanishes in the normalization employed here, see Eq. ([I]). The approximate 
formulas given above and in the following are all valid within 1% for 0.5 < rj < 0.6 and 
160 GeV < Mt < 170 GeV. Furthermore, the last line of the latter equation gives the 
numerical values of the individual contributions for rj = 0.56 and M t (M t ) = 165 GeV. All 
together, these 0(a) contributions lead to a shift of 0.00002 in CV(m&). Given that in 
our normalization CV(m&) ~ —0.0053 at NNLO QCD, this correction is quite small. The 
analogous expression for Cg(fi) reads 

Cg(/i) « -41.496 e-4-8790^0.4747 3.0.4738 + 75 53^-6.4866 ^0.493! ^.-1.9596 

+ 13.565 e - 6 - 0599 V- 3409 ( 15 ) 
« -4.0517 + 0.1579 + 0.1173. 

Here and in the following numerical results we have set fio = M w . We have distinguished 
between the three terms in Eq. ()13|) . Again, the last line in the above equation gives the 
numerical value of the three contributions for rj = 0.56 and M t (M t ) = 165 GeV: this elec- 
troweak correction to Cg(mb) is around —0.00226, thus much larger than the one in Eq. (jllj) . 
Moreover, it is dominated by the first term in Eq. (|13|). namely by the 0(a) matching con- 
ditions, and particularly by those of the electroweak penguin operators. Numerically, it 
corresponds to approximately —3% of the NNLO QCD coefficient. The corresponding 
expression for C w (fi) is 

CflM W -20.835 e -4.9135 ,^0.4957 x 0.8394 + g^g ^15.347^8.8810 

(16) 

w -3.3340 + 0.0720. 

The shift in Cio(rrib) is about —0.00195, and is again dominated by a term stemming from 
the higher order matching corrections. Like in the case of C-?([i), the term resulting from 
the leading order QCD running is exactly zero since there is no mixing of the four-quark 
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operators into Qio at 0(a s ), and the 0(a) correction to Cio(fio) is zero too. Numerically, 
it corresponds to approximately +3% of the NNLO QCD coefficient. 

The experience with radiative decays [9] shows that electroweak corrections to the 
matching conditions may end up being the dominant electroweak contributions to the low 
scale coefficients. In the notation adopted here, such contributions to CV(yUo), Cg(/io) and 
Cio(/-*o) are of 0(aa s ) and therefore formally enter at NNLO or 0(aa™L n ) in the b — > s£ + £~ 
amplitude. Although a complete analysis at this order is clearly beyond the scope of our 
paper, it is worth investigating the magnitude of the potentially dominant terms, such as 
electroweak effects in the matching enhanced by factors 1/ sin 2 9 W m 4.3 or by powers of 
M t . The general expression for C\ (/i) at the low scale is 

+ ^Wo^Guo) + 77tfj 1 Wo)<1 1) Guo) + u?\^^)cf\^) . 1 

The second, third, and fourth term in this equation are known completely, from [3] and 
the above discussion. The first term is partially known: in particular C^(hq) and Cg 2 g(/io) 
are known from [9], while only the part of Cfl(^o) and C^l (fio) proportional to Mf is 
known from [23]. We also include in the following the leading contributions to C^i (a*o) _ 

C^g (jUo) calculated in [3, 21] 5 . The last term in Eq. (JT7J) would require the knowledge 
of the three- loop 0(aa 2 s ) ADM and will be neglected in what follows. Again, rather than 
giving analytic expressions for these NNLO corrections, we will only give some approximate 
formulas and refer to a forthcoming publication [22] for more details. We use the results 
of [3,9,21,23] for the 0(aa s ) initial conditions of the relevant Wilson coefficients and find 
for the correction to CV(/i) 

Ci%) « 2.9457 e -0.4323^1.3783 x 0.9498 + l 2m ^5.4241, ^2.1365 ^-0.5479 

- 6.8887 x lo^e" 34 - 605 V 3 - 470 ^ 677 + 50.5 57 e' 5 - 9375 " r/ 2 ' 6489 x," ' 0736 (18) 
« 4.0718 + 0.0077 - 0.4315 + 0.3522 . 
The analogous expression for Cg(/x) reads 

eg (jj) « 2.3684 x 10 11 e -3i-006^i2.546 ^-o.gesi _ m52 ^mrir, ^2.3626 x -i.7mi 

+ 519.80 e - 6 - 5763 ^ 3 - 1984 x?- 4666 + 94.434 e" 4 - 8623 " r/ L1723 x° t 5775 (19) 
« 1.1803 - 0.3190 + 4.0021 + 7.2083 . 
Finally, for the 0(aa s ) correction to Cio(aO we § e t 

C%{p) « 0.3106 x 2 - 6056 - 242.56 e- 4 - 7109 ^ L5336 xr a7682 

+ 80.747 e -4.9425 ,^1.2373 x 0.7343 ( 20 ) 

g 13.144- 2.3657 + 7.1098 . 

5 Since the calculation of Ref. [21] has been performed in a different operator basis one has to perform 
a non-trivial change of scheme. See [22] for details. 
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In the last equation the contribution of the third term in Eq. ()17|) is exactly zero. In 
Eqs. ()18|) . (Jinj) and (f2Uj) we have set M H = 115 GeV and the last line always corresponds to 
the numerical value of the four individual contributions to Eq. ()17j) . employing rj = 0.56, 
M t (M t ) = 165 GeV, and setting /i = M w and /if = M t . For what concerns CV(mb), these 
0(aa s ) corrections are twice as large as the 0(a) correction, and lead to a shift of 0.00004. 
On the other hand, in the case of Cg(mb) and Cio(?rib) they are much smaller than the 
corrections of 0(a), as they amount to only 0.00016 and 0.00018. We stress that they are 
incomplete and we study them to check the potential size of the higher order contributions. 
If we add these partial higher order results to the complete leading and next-to-leading 
electroweak effects considered above, the values of the Wilson coefficients of Qj, Q$ and 
Q w are changed by —1.1%, —2.9% and +2.3%, respectively. 

At this point we also note that the electroweak corrections to Cio(/-tf>) depend rather 
strongly on the t-quark mass renormalization scale fi t . Using M t (M w ) = 175 GeV instead 
of M t (M t ) = 165 GeV, which usually leads to larger higher order QCD corrections, one 
finds a correction of +5.3%. This scheme dependence stems from the parts of C^^/io) 
proportional to Mf computed in [23], which we have evaluated in the MS scheme for the 
t-quark mass. On the other hand, the /^-dependence of the QED corrections to C 7 (/Xb) 
and C 9 (fib) is very weak. 

Let us now verify that the RGE has indeed preserved the correct normalization of the 
electroweak coupling for the matching contributions to Q 9 and Q w , as anticipated at the 
beginning of this section. To this end, we first consider the effect of the (3^ terms in the 
ADM of Eq. (JHHD- Using the explicit form of U^\/i, /io), one sees that they modify the 
contribution of Cg(/io) to the low-energy Wilson coefficients as follows 



U s%(H, ^o)C 9 (/i ) 



a 



a(/j,) 



U s °gg(^, /^o)C 9 (/i ) , 



Us%(fi, /io)C 9 (/i ) 

(21) 



where the equality holds up to small higher order terms in a and up to subleading loga- 
rithms in a s . Indeed, the consideration of subleading QCD effects, like those contained in 
U^\lJ,, /io), shows that QCD does not affect the proper normalization of the electroweak 
coupling, as anticipated above. Equation (|2T|) and the analogous relation for Cio(/io) 
demonstrate that the physically motivated normalization of the electroweak coupling is 
preserved by the RG evolution even for our choice of normalization of the operators: since 
the matrix element is multiplied by a(/i), the rescaling factor in front of Cg(no) brings the 
normalization of the electromagnetic coupling back to the electroweak scale. All other large 
logarithms described by the QED-QCD evolution matrices Ue l \fi, yUo) are solely induced 
by the running of the operators. 

While Us°\fi, /io) and (jf^ (/i, /io) are scheme- independent, Us^ (/i, /io) and U^ifa }J>o) 
depend on the renormalization scheme in a way that is detailed for example in Ref. [21]. 
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The scheme dependence is canceled by the 0(a s ) and 0(a) matrix elements of the effective 
operators. In fact, any calculation that goes beyond leading logarithms requires the knowl- 
edge of the one-loop matrix elements. However, we have seen above that the dominant 
effects originate at the weak scale and have often to do with 0(a/ sin 2 6 W ) contributions to 
the matching conditions. These contributions are all scheme- independent. One could also 
absorb some contributions to the QED matrix elements by replacing with a(m^), 
as discussed at the beginning of this section. However, numerically this would lead to a 
small difference and we will avoid it. We therefore expect our incomplete calculation to 
capture the dominant electroweak effects. Moreover, any complete calculation of the QED 
matrix elements should take into account real photon radiation, which is very sensitive to 
the details of the experimental set-up and requires a dedicated study at the experimental 
level. 



4 Phenomenology 



We are now ready to combine all our results and update the prediction for BR«. We first 
write the integrated branching ratio as in [3,4]: 



/•0.25 i 

BR« = BR[B - XM / ds _ — 



dY[B — > X s £ + £~] 



ds 



(22) 



where 



dT[B -> X s £+£~] _ Glml polc \V*V tb f ( a(m b )\ 2 _ 2 



ds 



+ (1 + 25)1 & 9 % R (s) 



48vr 



An 




(23) 



and the effective Wilson coefficients include all real and virtual QCD corrections [4] as in- 
dicated by the index BR, whereas the last term encodes finite bremsstrahlungs corrections 
computed in the second paper of Ref. [4]. At 0(a) the definition of the effective Wilson 
coefficients is affected by the presence of the electroweak penguin operators Q^-Qq and 
of Q\ — for a detailed discussion we refer to [22]. In order to cancel the strong depen- 
dence on the bottom quark mass in the factor rn|, it is customary to normalize BR^ to 
the experimental value of the Branching Ratio (BR) for the inclusive semileptonic decay 
BR[i? — > X c tv\. However, this normalization introduces a strong dependence on the charm 
quark mass, which is not known very accurately. The ensuing uncertainty is about 8% [5]. 
An alternative procedure [25,26] consists in normalizing the B — > X s £ + £~ decay width to 
T[B —>■ X u £u], and then to express BR[i? — > X u £v\ in terms of BR[B — ► X c £u] and of the 
ratio 



C 



ub 



do 



T[B -»■ X c £v\ 
T[B X u £v\ 



(24) 
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S S 

Figure 5: Scale dependence of the differential BR« using the different normalizations of 
Eqs. (I22|) and (I25|) in the partial NNLO QCD approximation (Ihs). The discontinuity at the 
end of the low-s region reflects the fact that the 0{a 2 s ) matrix elements of Q\ and Q 2 have 
not been implemented for s > 0.25 [8], while the divergence in the J/ip region is due to the 
1/m 2 corrections. Furthermore the relative size of the NNLO and electroweak (EW) effects 
introduced in this paper are shown when sequentially included (rhs). 

which can be computed with better accuracy [25]. In other words, we will use 

V ub 



BR[B -> XM 
otiee — ^ 



V, 



cb 



1 dv\B^x s i + n 

o. 05 T[B^XJv] ds 1 ; 



The main uncertainty in C comes from non-perturbative effects. Using an updated deter- 
mination of the heavy quark effective theory parameter Ai = (—0.31 ±0.10) GeV 2 [27], and 
following otherwise Ref. [25], we obtain C = 0.581 ± 0.017. For the inclusive semilep- 
tonic b -> c branching ratio we use BR s/ = BR[B -> XJu] = (10.74 ± 0.24)% ob- 
tained from the latest Heavy Flavor Averaging Group results [28] after subtracting the 
b — > u transition. For the Cabbibo-Kobayashi-Maskawa (CKM) factor \V t * V t b\ 2 /\ V c b\ 2 
we use 0.969 ± 0.005 [27], but we actually perform the calculation keeping top, charm 
and up quark sectors separated [3]. Finally, we employ the c- and 6-quark pole masses, 
m CjPO i c = (1.4 ± 0.2) GeV and m feiPolc = (4.80 ± 0.15) GeV. 

The use of the 6-quark pole mass in Eq. (|23|) and in the calculation of the semileptonic 
width leads to large perturbative QCD corrections both in the numerator and denominator 
of Eqs. (|22j) and (J25|) . Since this is an artifact of the choice of scheme for the bottom quark 
mass, they clearly tend to cancel in the ratio. As a second improvement with respect to 
previous analyses, we keep terms through 0(a 2 ) in the denominator and expand the ratio 
in Eq. (|25j) in powers of a s . By making explicit the cancellation of large contributions, 
we improve the convergence and stability of the perturbative series, as we have explicitly 
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partial NNLO 


tt(2) 
(J s,92 


a (2) 

u 99 


LO + NLO EW 


EW 


BR« as in 


1 c;nq +0 - 035 
J..JUy_o,o88 


1 coq+0.030 




l. / ±UO_ .087 


1 449+ ' 054 


BR« as in <|27j) 


1 64fi +0 - 038 


1.668lo;o8? 


1 ^07+0-041 
l.OJ 1 -0.087 


1 cqo+0.060 
l.dOZ_ 094 


1 c,7_+0.059 
i -°' 4 -0.075 



Table 2: BR« integrated over the low-s and low-g 2 region in units 10 6 . The different correc- 
tions are added in sequence. 

checked. Like in previous analyses, due to the peculiarity of the perturbative expansion 
for b — > s£ + £~, we retain some large and scheme-independent higher order term in the 
amplitude squared: for instance the terms |C*7 ff (s) | 2 and \C^ (s)| 2 in Eq. are expanded 
up to 0(a s ). In addition all known 0(a s ) terms of \Cf{s)\ 2 and Re{Cf{s)Cf{s)*) are 
kept. Our procedure differs slightly from the one proposed in the first paper of [4], which 
consists in absorbing the numerically small coefficient Cg° g (fib) into C^Q^), and in sub- 
sequently performing the a s expansion of the numerator dT[B — > X s £ + £~]/ds only. In the 
partial NNLO approximation, that is, setting t/ s 92 (/X6, Hq) and a 99 (s) to zero, this alterna- 
tive method yields somewhat lower values for the normalized differential rate, well within 
the scale dependence of our results. 

We also include power corrections of 0(l/m 2 ) [29] and 0(l/m 2 ) [30] to the differential 
decay rate of B — > X s £ + £~ as well as 0(l/mf) corrections to T[B — > X c £u] and T[B — > 
X u £v\ — see for example [31] - - and expand the normalized differential decay rates given 
in Eqs. (122)1 and (|25)l in inverse powers of nib and m c . 

The partial NNLO results before inclusion of the two contributions discussed in Section 
2 are shown in the left plot of Figure |_I where the bands include only the residual scale 
dependence, calculated performing a scan in the ranges: [40-120] GeV for the matching 
scale in the charm and up sector, [60-180] GeV for the matching scale in the top sector, 
[2.5-10] GeV for the low-energy renormalization scale. In the last case, the central value 
is set to 5 GeV. As can be nicely seen by comparing the gray with the yellow band, the 
residual scale dependence for low s is about 50% smaller in our approach than without 
the two improvements introduced between Eqs. ()23|) and (|23j) . On the other hand, the 
scale uncertainty is practically unaffected by the new NNLO contributions — not shown 
in Figure El 

The results for the integrated BR are given in the first line of Table |2l where the 
first column refers to the partial NNLO calculation according to [4] but implemented as 
described above, the second column to the same calculation supplemented by our result 
for U[ ; 92 (/ife, /io) ; the third column to the complete NNLO, that is, including the two- 
loop QCD matrix element a^ 9 \s). The fourth column refers to the NNLO calculation 
including LO and NLO electroweak effects, and the fifth column to the QCD computation 
supplemented by the electroweak effects through NNLO calculated in the previous section. 
Furthermore, in the last two columns the dominant QED corrections to the semileptonic 
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decay amplitude [32] have been included. If we had only rescaled the matching conditions 
for Q Q and Q w and for the penguin operators by a factor a(M w ) / ot{m h ) as discussed in the 
previous section, the shift would be close to +7%. This effect is overcompensated by other 
electroweak corrections of around —7% and by the QED correction to the semileptonic 
rate, which by itself amounts to around —1.5%. We would like to stress that this is an 
unexpected accidental cancellation. The right plot of Figure El shows A[dBR u (s) / ds] 
the differential BR normalized to the partial NNLO result — where we have sequentially 
included the corrections listed above. 

Including the main parametric uncertainties, our final result reads 

BR U (0.05 < s < 0.25) = (26) 

[1.442 ±°;™» \ Mt ±°;°*| | scale ± .041 c ± 0.032 BRs; ±° Z \ mc ± 0.002 m J x 10~ 6 . 

Here the whole dependence on Ai is absorbed into the factor C. As already in the case of 
C, we neglect effects suppressed by three powers of the heavy quark masses in B — > X s £ + £~ 
[33], which in the low-s region should not exceed 2% anyway. The uncertainty related to 
the CKM matrix elements is below 1%. To be conservative, one could also consider an 
additional 2% uncertainty due to unknown electroweak effects. Combining in quadrature 
all uncertainties given in Eq. ()26|) with the latter three, the total error is about 9%, very 
similar to the one for inclusive radiative decays [25]. 

Though useful for comparison with previous analyses [4], the BR integrated in the region 
0.05 < s < 0.25 is an idealization. A quantity closer to experiment is the BR integrated 
in the region 1 GeV 2 < q 2 < 6 GeV 2 [5] . The results for the integrated BR employing 
this range are given in the second line of Table 121 together with their scale dependence. 
Unfortunately, the choice of the integration variable introduces an enhanced dependence 
on the bottom quark mass. In our calculation we have employed the 6-quark pole mass, 
which is subject to a much larger uncertainty than an appropriate short- distance 6-quark 
mass. For example the IS or the kinetic low-energy 6-quark mass [27] are known to better 
than 50 MeV. Our final result is 

BR U (1 GeV 2 < q 2 < 6 GeV 2 ) = (27) 

[1.574 ±oioo |m« ±ao67 \m b ±ao75 Iscaic ± 0.045 c ± 0.035 B R a! ±aoi3 Lj x 10 6 . 

The total error is about 10% and is again dominated by the uncertainty on the top quark 
mass, which will soon be reduced by a factor of two by CDF and DO at the Tevatron. 
Moreover, the substantial error from the 6-quark mass is an artifact of the employed scheme, 
which can be reduced by a factor of three or even more by changing the renormalization 
scheme of the fe-quark mass. Hence it should not be interpreted as a limitation. 

For what concerns the FB asymmetry, we have investigated the impact of the choice of 
normalization and of the electroweak corrections on the so-called unnormalized asymmetry 
defined as in [4] 

„ BR[B^XM f 1 , „ d 2 Y[B -> X s e + £-\ , „ , 

A FB (s) = — I dcos9 e l —— '- isgn(cos^), (28) 

1 [B — > X c lv\ J_i ds d cos Ve 
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so without 0(l/m^ c ) 


s with 0(1/ ml c ) 


unexpanded A FB (s) as in (J2HI) 


n 1 fi n+0.008 
u.iuu„ 005 


n 1 R7+0.008 
u - io ' -0.006 


expanded A FB (s) as in (|3L)j) 


n ic- 1+0.005 
U.lOl„ 004 


n i c; fi +0.005 

U.lOU_ 005 


expanded A FB (s) as in (|31j) 


n 1 c;7+0.005 
u - iJ ' -0.009 


f) 1 fi Q+0.007 
U.lDO_o,010 


unexpanded A FB (s) as in (J2SJ) + EW 


n i fi4+o.oo4 

U.lU / ±_ 004 


n 1 7n+o-°o 4 

U - L ' u -0.005 


expanded A FB (s) as in fl3Ul) + EW 


n i p-4+0.004 
U.lO^_ 004 


n i c;q+0.004 

U.ldJ_ 004 


expanded A FB (s) as in (pTTf + EW 


n i fi Q+0.006 
U.lOU_ 009 


n ififi +0 - 007 

u.ioo_ oio 



Table 3: Zero of FB asymmetry before and after the inclusion of \jm\ and l/m 2 c power 
corrections. 



where 
-i 



i 



rf cos 6> £ — sgn(cos 6^ 

as a cos 6^ 



G fM m b, P oie \ v ts v tb[ I a{m b 



48vr 3 



47T 



(29) 



6 Re (CSbC^CSjb^)*) - 3s Re (Q% B (s)C^ FB (sy ) + 



and the effective Wilson coefficients include all real and virtual corrections calculated in [4] 
as indicated by the subscript FB. The last term of Eq. (|29|) encodes finite bremsstrahlungs 
corrections evaluated in the last reference of [4] which we have not included 6 . At O(a) the 
definition of the effective Wilson coefficients is affected by the presence of the electroweak 
penguin operators Q^-Qq and of Q\ — for a detailed discussion we again refer to [22]. 

Implementing all the relevant NNLO QCD corrections we obtain for the position of 
the zero of the unnormalized FB asymmetry (|2*8*|) the values given in the first column of 
the first line of Table El in agreement with [4,8]. While in Refs. [4,8] only the low-energy 
renormalization scale was varied, the errors in Table El also include the matching scale 
dependence. 

Since the position of the zero of the FB asymmetry is known to be especially sensitive 
to physics beyond the SM, a comprehensive study of the residual theoretical error attached 
to it is important. In order to get a grasp for the size of higher order corrections, we 
consider the normalization to the charmless semileptonic decay rate 



.4 



fb(s 



BR[5 -> X c £u] 1 
T[B^XJv] C 



ab 



V, 



cb 



a cos Ui - 

i 



ds d cos 9 



- sgn(cost^ 



(30) 



in combination with C and perform an expansion of the ratio of Eq. ()30|) in a s to improve 
the convergence of the perturbative series. The corresponding result is displayed in the 



3 These corrections turn out to be below 1%. 
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Figure 6: Unnormalized FB asymmetry A FB (s) in the low-s region including l/mf and l/m 2 c 
power corrections. Compared is the NNLO QCD result for the normalization given in Eq. ()28|) 
and the one given in Eq. (I30|) performing a subsequent expansion in a s . Furthermore we also 
display the result of Eq. (J3UJ) after the inclusion of electroweak effects (central yellow band). 



first entry of the second line of Table El Note that in this case the numerical value obtained 
for So does no longer correspond to the zero of the curly bracket in Eq. (|29|). The difference 
between the values of §q in the first and second line of Table 01 is obviously due to different 
higher order terms and can give an indication of the residual theoretical error. 

We also consider the so-called normalized FB asymmetry, which is certainly a quantity 
closer to experiment: 



A FB (s) 



dT[B — > X s £ + £~]/ds 



a cos &£ - 

i 



Xl+t 



ds d cos 9/ 



sgn(cos 9i 



(31) 



Given the small value of the LO term of dT[B — > X s £ + £~]/ds, a naive a s expansion of 
this ratio is problematic. However, the a s expansion converges reasonably well if the 
perturbative expansion is reorganized following the first paper of [4], namely absorbing 



Cfg(fib) into C^g (/i;,). With respect to the unexpanded form, which obviously reproduces 
the value given in the first column of the first line of Table El the zero is shifted toward 
lower values, as can be seen by comparing the first and the third entry in the first column 
of that table. Finally, the next three entries of the first column of Table El correspond to the 
results for §q in the three approaches described above after the inclusion of the electroweak 
effects. 
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Figure 7: Normalized FB asymmetry A FB (s) in the low-s region in the NNLO QCD approxima- 
tion as given in Eq. (13 ip without expansion in a s . Furthermore the effects of including power 
corrections (central band) and higher order electroweak effects (lower yellow band) are shown. 

Unlike previous analyses, we also include power corrections of 0(1/ mf) [29] and 0(1/ m 2 c ) 
[30] in the computation of the zero of the FB asymmetry. In doing so we always expand 
the asymmetry in powers of the heavy quark masses. The corresponding numerical values 
for s can be found in the second column of Table H3 Comparing the results in that table, 
one observes that higher order electroweak effects amount to around +2%, independent on 
which normalization is used and on whether an expansion in a s is performed or not. The 
0(1 /ml) and 0(l/m 2 c ) corrections turn out to be at the level of +4%. The unnormalized 
FB asymmetry Afb(s) calculated in the two above approaches in the vicinity of the zero 
is shown in Figure |H1 whereas in Figure [7| we plot the unexpanded Afb(s), showing the 
importance of the power and electroweak corrections. 

A glance at Table |3] and Figures |H1 and [7| shows that the dependence of the zero on 
the specific method used to calculate it is generally larger than the scale dependence 
associated with it. This suggests that the observed scale dependence alone does not provide 
a reasonable estimate of higher order QCD and electroweak corrections, as pointed out 
already in the third paper of [4]. Following a reasoning similar as for BR«, we adopt as 
the central value for §o the average of the two expanded results in Table H3 and assign to it 
a rather conservative total error of 6%, that covers the whole range of possible values: 



The uncertainties due to M t and the non-perturbative parameters Ai and A2 are at the 0.1% 



s = 0.163 ±0.010 theory . 



(32) 
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level or less. Converting this result into a value of q 2 introduces again a large uncertainty 
due to the bottom quark mass. We find 

q 2 = (3.76 ± 0.22 thcory ± 0.24 m J GeV 2 , (33) 

with a total error close to 9%. However, also in this case the error due to the 6-quark mass 
can be drastically reduced by performing the calculation in a different mass scheme. 

5 Summary 

Rare semileptonic B decays are going to play an important role in the search for new physics 
at the B factories and upcoming flavor physics experiments at the Tevatron and the LHC. 
From the theoretical point of view the inclusive mode at low q 2 is particular interesting, 
since it can be accurately computed in the SM. In this paper we have first completed 
the NNLO QCD analysis of B — > X s £ + £~ with the calculation of the last two missing 
components, and then discussed electroweak effects. The new NNLO QCD contributions 
are relatively small — they amount to around +1% and —4% — and tend to partially 
cancel each other. Our treatment of electroweak effects is complete only for what concerns 
the leading logarithms of photonic origin but, pending an exhaustive investigation of the 
missing matrix elements [22,24], we expect it to capture all the dominant contributions. 
The impact of electroweak effects depends on the scale employed for a(n) in the prefactor 
for BRa. As a result of accidental cancellations, they change the value of the low-energy 
BR by only —1.5% for fi = 0(mf,), while if one adopts /i = 0(M W ) the effect is much more 
relevant, about —8.5%. Our calculation drastically reduces the large uncertainty of almost 
8% related to higher order electroweak effects, which has affected previous analyses. We 
have also calculated the FB asymmetry. The location of its zero is shifted by approximately 
+2% due to the electroweak corrections. 

We have also introduced a few improvements in the calculation of the decay rate in the 
region of low dilepton invariant masses, aimed at reducing the theoretical error. After a 
careful study of the various sources of uncertainty, we estimate a 9% to 10% error on our 
result for the BR in the s = 0.05-0.25 or q 2 = 1-6 GeV 2 window, the dominant source 
of which is the uncertainty on the top quark mass. For what concerns the zero of the FB 
asymmetry, our estimate for the residual uncertainty is 6% to 9%. 
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Appendix 

We report below the complete 0(a s ), 0{a 2 s ), 0(a) and 0(aa s ) ADMs making explicit 
their dependence on the coefficients of the QED and QCD (3 function: 
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(34) 



For five active quark and three active lepton flavors the values of the MS 8 coefficients are 
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Adopting the ordering of operators as introduced in Eq. (JT|) the matrices 7! and 7J 
take the following form 
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The one-loop 0(a s ) mixing of the electroweak penguin operators into Qg, as well 

as the mixing of Q\ has never been computed before. The same applies to the complete 
two- loop 0(a 2 s ) mixing 7 of Q®-Q® and Q\. As far as the remaining entries are concerned 
our results agree with the well-established results in the literature [3,10,12,34,36]. 

The matrices 7e°' ) and 7! are given by 
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7 The self-mixing of Q^-Q® has been computed in the so-called "standard" operator basis in [35]. 
Since the operator basis of the latter papers differs from the one used here one has to perform a non-trivial 
transformation. This can serve as a cross-check of the direct calculation. See Ref. [22] for details. 
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In only the entries corresponding to the self-mixing of Qg and Q\o are new, while all 
other elements agree with those of [10]. On the other hand, -yi is entirely new. 

The two-loop ADMs 7! and 7! depend on the definition of evanescent operators. We 
have followed the definition employed in [12,34] extending it to include all electroweak and 
semileptonic operators: 

Ei = (s L l,lvl P T a c L )(c L YYYT a b L ) - 16Q! , 

E 2 = {sLl^u1 P c L ){c L Yl u Yb L ) - I6Q2 , 

£3 = (sLl.lvlpl.lrb^E^YYY^q) + 64Q 3 - 20Q5 , 

£4 = (s L ^ lulplalr T% L )J2 q (qYYYri T T a q) + 64Q 4 - 20Q 6 , 

£ 3 Q = (sLl»lul P l«lrb L m q Q q ((lYl u l P l a l T q) + 64Q? - 20Q« , 

E? = (sL7M7,7p7 CT 7rT°& i )EA(?VV7V7 r 7 10 ?) + 64Q 4 Q - 20Q« , 

(40) 

E\ = ±(s Ll ^ ulplalT T a b L )(bYYl P ri T T a b) + f(s LllI T a b L )(h^T a b) 
- l{s Ll , lvlp T% L ){bn u YT a b) - 2Q\ , 

E\ = (s L l»lvl P l<rlrb L )(W L l v Yl a l T b) - f {SLl P b L )(W L b) 
+ l(s L ^ ulp b L ){bYYYb) - 256Q? , 

e 9 = ^(s Llplulp b L )j: e (lYYY^ - 10Q9 + 6Q 10 , 
£10 = $(s Llplulp b L )j: e (iYYiV) + 6Q 9 - 10Q10 • 

At this point a comment concerning the two-loop mixing of Q\ is in order. Using the 
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definition of Q\, E\ and E\ as given in Eqs. (JQ and (|4U|). the electromagnetic operator 
Q\ mixes into new physical operators beyond the one-loop level. Needless to say, such a 
mixing does not alter the results presented throughout the paper. Further details on the 
two- loop mixing of Q\ will be given in a forthcoming publication [22]. 
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